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We present a one-loop analysis of the pairwise production of neutral Higgs bosons (h A ,H°A ) 
\ at linear colliders, such as the ILC and CLIC, within the general Two-Higgs-Doublet Model 



(2HDM). We single out sizable radiative corrections, which can well reach the level of |5 o\/o ~ 
50% and may be either positive (typically for y/s ~ 0.5 TeV) and negative (for s/s > 1 TeV). These 
large quantum effects, obtained in full agreement with the current phenomenological bounds and 
the stringent theoretical constraints on the parameter space of the model, can be traced back to 
the enhancement capabilities of the triple-Higgs self-interactions - a trademark feature of the 
2HDM, with no counterpart in e.g. the Minimal Supersymmetric Standard Model. In the most 
favorable scenarios, the Higgs-pair cross sections may be boosted up to barely 30 fb at the fiducial 
center-of-mass energy of 500 GeV - amounting to 10 3 ) events per 500 fb~' of integrated 
luminosity. We also compare these results with several complementary double and triple Higgs- 
boson production mechanisms at order ff(a^ w ) and leading ^(a^ v ), and we spotlight a plethora 
of potentially distinctive signatures of a Two-Higgs-Doublet structure of non-supersymmetric 
nature. 
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1. Introduction 

The Two-Higgs-Doublet Model (2HDM) is a particularly simple extension of the Standard 
Model (SM) which already encompasses outstanding new phenomenology [1]. In addition, a Two- 
Higgs-Doublet structure naturally emerges as the low-energy realization of some more fundamental 
theories - viz. the Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) [2]. 
The 2HDM spectrum contains two neutral CP-even (h°, H°), one CP-odd (A ) and two charged 
Higgs bosons (H*) and, as a most distinctive feature, it allows triple (3H) and quartic (4H) Higgs 
self-interactions to be largely enhanced - in contrast to the MSSM, wherein such couplings are 
restrained by the gauge symmetry. The phenomenological impact of such potentially large 3H 
self-interactions has been actively investigated at linear colliders within a manifold of processes, 
to wit: the tree-level production of triple Higgs-boson final states [3]; the double Higgs-strahlung 
channels hhZ° [4]; and the inclusive Higgs-pair production via gauge-boson fusion [5]. Likewise, 
the 77 mode of linac facilities has also been considered, in particular within the loop-induced 
single [6] and double Higgs production processes [7]. All these mechanisms would be capable to 
yield large Higgs production rates and furnish experimental signatures which, owing to the clean 
environment inherent to linear colliders, might enable for precise measurements of e.g. the Higgs 
boson masses, their couplings to fermions and gauge bosons and their own self-interactions; this 
means, to reconstruct the Higgs potential itself. Moreover, these large rates could be also revealing 
by themselves, as they could not be deemed e.g. to a SUSY origin due to the intrinsically different 
nature of Higgs self-interactions. 

Similary, double Higgs boson (2H) production may be also instrumental at future linac facili- 
ties [8]. Such processes cannot proceed at the tree-level in the SM, which means that, if we would 
detect a sizable rate of 2H final states e.g. of the sort e + e~ — > h°A ;H () A°;H + H~, this would entail 
an unmistakable sign of new physics. Nonetheless, a tree-level analysis of these processes is most 
likely insufficient to disentangle e.g. SUSY and non-SUSY extensions of the Higgs sector, as the 
tree-level &A°Z couplings (h = h°,H°) are purely gauge, and hence both models can give rise to 
similar cross-sections at the leading-order. One-loop corrections to these 2H production rates are 
therefore required to this effect. A number of studies are available within the MSSM [9] whilst, on 
the 2HDM side, the efforts were first concentrated on the production of charged Higgs pairs [10] 
and, only very recently, they have been extended to the neutral sector [11]. 

2. Computation setup 

The general 2HDM [1] is obtained upon canonical extension of the SM Higgs sector with a 
second SUl{2) doublet with weak hypercharge Y = +1. The seven free parameters A,- in the general, 
CP-conserving, 2HDM can be sorted out as follows: the masses of the physical Higgs particles 
(Mffi, M H o, M A o, M H ±), tan/3 (the ratio of the two VEV's (H?)), the mixing angle a between the 
two CP-even states; and one Higgs self-coupling A5. Additionally, the Higgs couplings to fermions 
must be engineered so that no tree-level flavor changing neutral currents (FCNC) are allowed: this 
gives rise to the following basic scenarios, namely Types-I and II 2HDM - see [1] for details. 
Further constraints must be imposed to assess that the SM behavior is sufficiently well reproduced 
up to the energies explored so far, namely: i) the perturbativity and unitarity bounds [12]; ii) 
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the approximate SU(2) custodial symmetry, which requires \8p2hdm\ ^ 10 3 [13]; and Hi) the 
consistency with the low-energy radiative B-meson decays (which demands M H ± > 300 GeV for 
tan/3 > 1 in the case of type-II 2HDM [14]). We refer the reader to [11] for further details on the 
model setup and the constraints. 

Hereafter, our endeavor will be basically threefold: i) to single out the regions in the 2HDM 
parameter space for which the e + e~ — > A h°/A () H () production is optimized; ii) to quantify the 
numerical impact of the associated quantum effects; and iii) to correlate the latter with the en- 
hancement capabilities of 3H self-interactions. For denniteness, let us concentrate on the h°A° 
channel. Our starting point shall be the e + e~ — > h°A° scattering amplitude at 1-loop, 

where \J Z h o h o stands for the finite WF renormalization of the external h° leg, which we shall ex- 
pand as yj Z h o h o = 1 — jReL' h o hQ (M^ ), in order to retain only those contributions which are at 

leading-order in the triple-Higgs self-couplings. In turn, + -^A°h° i ncuj des the entire set of 
ff(a^ w ,a ew a em ) one-loop contributions, to wit: i) the vacuum polarization corrections to the Z°- 
boson propagator and the Z° — y mixing; ii) the loop-induced yA°h° interaction; iii) the vertex 
corrections to the e + e~Z° and h°A°Z° interactions; and iv) box-type contributions. Finally, the 
1-loop counterterm S^l_ AOh0 guarantees that the overall amplitude M is UV-finite. Such 
counterterm ought to be anchored by a set of renormalization conditions; for the SM fields and 
parameters, we stick to the conventional on-shell scheme in the Feynman gauge, see e.g. [15]. 
Concerning the renormalization of the 2HDM Higgs sector, we also employ an on-shell scheme 
whose implementation is discussed in full detail in [1 1]. 

3. Numerical results 

The basic quantities of interest are: i) the predicted cross section at the Born-level cy( ) and at 
1-loop cj( 0+1 \ in which we include the full set of G{d^) corrections, and also the leading G(a% v ) 
which come from the squared of the scattering amplitude and ii) the relative size of the 1- 
loop radiative corrections, which we track through the parameter 8 r = g( ° + ^ {T( ° > . Throughout the 
present work we make use of the standard computational packages [16]. For definiteness, we sort 
out the Higgs boson masses as follows: 

M h o[GeV] M H o[GeV] M A o [GeV] M H ± [GeV] 
Set A 130 150 200 160 

SetB 150 200 260 300 

We shall use type-I 2HDM for Set A, and type-II for Set B - as, in the latter case, M H ± is suitable 
to elude 3§{b — > sf) constraints. The predictions for both sets are quoted in Table [I] for different 
values of the /jA°Z° tree-level coupling and maximally enhanced 3H self interactions - namely 
tan/3 = 1 and X 5 ~ -10 (resp. -8) for Set A (resp. Set B). 

These results single out large quantum effects (\Sa\/a ~ 20 — 60%), which can be either 
positive (for y/s ~ 0.5 TeV) or negative ( > 1 TeV) and turn out to be fairly independent on the 
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h°A° 


H°A° 


^ = 0.5 TeV 


a = j8 a = — 7t/6 a — n/2 


a = P-7t/2 a = j3-7i/3 a = 


Set A 


Gmax [fb] 

S r [%] 


26.71 20.05 13.10 
31.32 31.42 28.81 


26.32 17.53 10.29 
48.45 31.82 16.10 


SetB 


Gmax [ft>] 
Sr[%] 


11.63 9.08 6.36 
35.17 40.68 47.86 


5.00 3.28 1.95 
71.81 50.22 34.15 




■s/s = 1 TeV 


a = j8 a — j8 — k/6 a — n/2 


a = P-n/2 a = p-7i/3 a = 


Set A 


Omax [fb] 
Sr[%] 


4.08 2.70 1.56 
-58.42 -63.28 -68.11 


3.59 2.52 1.55 
-62.62 -65.09 -67.75 


SetB 


Omax [fb] 
S r [%) 


6.11 4.22 2.86 
-30.16 -35.62 -34.58 


5.17 3.98 2.75 
-36.80 -35.14 -32.80 



Table 1: Maximum total cross section o"( 0+1 ) (e + e~ — > A°h°) and relative radiative correction (S r ) at y/S = 0.5 TeV and 
y 7 ? = 1 TeV, for Sets A and B of Higgs masses. The results are obtained at fixed tan ft = 1 and different values of a, 
with X$ ~ -10 (resp. -8) for Set A (resp. Set B). 



Set A 




12345 12345 
tan tan /3 



Figure 1: Contour plots of the quantum corrections S r (in %) to the e + e~ — ¥ A°h° cross-section as a function of tan/3 
and A5, for Set A of Higgs boson masses, a = j5 and ^/j = 0.5 TeV (left panel), *Js = 1 .OTeV (right panel). The shaded 
areas are excluded by the vacuum stability (top purple area) and the unitarity bounds (bottom yellow area). 



details of the Higgs mass spectrum, the particular type of 2HDM and the specific channel under 
analysis (h°A ,H°A ). The corresponding cross-sections at 1-loop lie in the approximate range 
of 2 — 30fb for yfs = 0.5 TeV - up to barely 10 3 — 10 4 events per 500fb _1 . As for the radiative 
corrections themselves, in Figure [l] we display their detailed behavior (for h°A°) as a function 
of tan/3 and A5, as well as their interplay with the unitarity bounds (lower area, in yellow) and the 
vacuum stability conditions (upper area, in purple). Notice that the former disallows simultaneously 
large values of tan/3 and X5, whereas the latter enforces A5 < 0. The largest attainable quantum 
effects (positive or negative, depending on ^/s) are identified in a valley-shaped region at tan /3 ~ 1 
and I A5 1 ~ 5 — 10. In such regimes, a subset of 3H self-couplings becomes substantially augmented 
- their strenght growing with ~ \X%\ - and stand as a preeminent source of radiative corrections - 
via Higgs-boson mediated 1-loop corrections to the /iA°Z° vertex. As a result the latter interaction, 
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which is purely gauge at the tree-level, becomes drastically promoted at the 1-loop order. Upon 
simple power counting arguments and educated guess we may barely estimate the loop-corrected 
h°A°Z° coupling as ~ T^X^/lSx 2 s x /(M^/s^/s), where f{M^/s,N^/s) is a 

dimensionless form factor. 




Figure 2: Total cross section cr(e + e~ — s> A°h°) (in fb) at the tree-level and at one-loop for different values of A5 and 
for Sets A (top panels) and B (bottom panels) of Higgs boson masses. Displayed is the behavior of o" as a function of 
^fs (left and center) and M h o (right panels). 



All these trademark phenomenological features become also patent in Figure |2], in which we 
explore the evolution of the h°A° production cross-section for the following setups: a-b) as a func- 
tion of y/s, for a = j5 (maximum h°A°Z coupling); c-d) the same, for a = n/2 (the so-called 
fermiophobic limit for the h° boson within type-I 2HDM); and e-f) as a function of M h o, for a = j8. 
The results are presented for Sets A (upper row) and B (lower row). The phenomenological foot- 
print of 3H self-couplings is easily read off from the dramatic differences in the one-loop produc- 
tion rates when the value of A5 is varied in the theoretically allowed range. As for the behavior 
with M h o, these plots illustrate a twofold effect on a; one is purely kinematical, and accounts for 
the decrease of the cross-section owing to the reduction of the available phase space; the other is 
dynamical, and explains why the tree-level cross-sections decouple faster than their corresponding 
one-loop counterparts; indeed, several 3H self-couplings get boosted with M h a and can partially 
counterbalance the phase space suppression. 
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0.5 TeV 


1.0 TeV 


1.5 TeV 


0.5 TeV 


1.0 TeV 


1.5 TeV 


Set A 


SetB 


CT(h U A°) [fb] 


26.71 


4.07 


1.27 


11.63 


6.11 


2.52 


CT(h°H°A ) [fb] 


0.02 


5.03 


3.55 


below thres. 


1.25 


1.33 


cr(H H+H-)[fb] 


0.17 


11.93 


8.39 


below thres. 


0.69 


2.14 


CT(h°h° + X) [fb] 


1.47 


17.36 


38.01 


0.92 


9.72 


23.40 



Table 2: Comparison of the predictions for the cross sections corresponding to several Higgs-pair and triple-Higgs 
production channels, for Sets A and B of Higgs masses; tan/3 = 1, a = j8, and three different values of the center-of- 
mass energy. The complementarity between the different channels at different energies becomes manifest here. 

4. Discussion and conclusions 

We have undertaken a comprehensive study of the neutral 2H production from a generic 
2HDM in the context of the ILC/CLIC colliders. The upshot of our analysis spotlights very sig- 
nificant radiative corrections, typically up to \8a/a\ ~ 50%, which may be either positive (for 
yfs ~ 0.5 TeV) and negative (for yfs ~ 1 TeV and above), and basically concentrated in the regions 
with tan/3 ~ 1 and |As| ~ 10, with A5 < 0. Such enlarged quantum effects can be basically traced 
back to the Higgs-mediated 1-loop corrections to the hfiPZ interaction, these being sensitive to po- 
tentially enhanced 3H self-couplings - a genuine feature of the 2HDM, with no counterpart in the 
MSSM. In the most favorable scenarios, the predicted 2H rates furnish a few dozen fb per 500 fb -1 
of integrated luminosity. In practice, these 2H processes would boil down to A — > bb/r + T~, and 
h° -»■ bb/r+T- or h° -»■ VV -»■ 41,21+ mising energy, depending on the actual value of M h o; all 
these signatures being feasible in the clean environment of linac facilities. Besides, the analysis 
of such pairwise Higgs events should be useful to enlighten the inner structure of the underlying 
Higgs sector, in particular if combined with complementary Higgs production channels such as 
e + e~ — > hhh [3] and e + e~ — > hh+X [5]. In Table || we quantify the latter statement by plugging 
(7(0+!) for h°A° together with the predicted (leading-order) production rates for several complemen- 
tary multi-Higgs channels in e + e~ collisions, assuming maximum 3H self-coupling enhancements 
(as e.g. Table |l|). This pattern of signatures at different y/s is highly distinctive of the 2HDM and 
could not be attributed to e.g. the MSSM, as they critically depend on the 3H self-couplings - 
which are fully inconspicuous in the latter model owing to supersymmetric invariance. Additional, 
and very valuable, information can also be obtained from the study of the quantum effects on the 
production of a neutral 2HDM Higgs boson in association with the Z boson [17]. 

We conclude that 3H self-couplings may stamp genuine footprints on multi-Higgs production 
processes, either at leading-order or through quantum corrections, and provide a plethora of com- 
plementary signatures whose detection should be perfectly feasible at future linac facilities and, if 
ever observed, might constitute a strong hint of non-standard, non-supersymmetric Higgs physics. 
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